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We reexamine the Glauber model and calculate the total reaction cross 
T— I ■ 

(<~^ ■ section as a function of energy in the low and intermediate energy range, where 

■ many of the corrections in the model, are effective. The most significant 



effect in this energy range is by the modification of the trajectory due to 
the Coulomb field. The modification in the trajectory due to nuclear field 
is also taken into account in a self consistent way. The energy ranges in 



2 ■ which particular corrections are effective, are quantified and it is found that 

when the center of mass energy of the system becomes 30 times the Coulomb 
j_j ■ barrier, none of the trajectory modification to the Glauber model is really 

required. The reaction cross sections for light and heavy systems, right from 
near coulomb barrier to intermediate energies have been calculated. The 
exact nuclear densities and free nucleon-nucleon (NN) cross sections have been 
used in the calculations. The center of mass correction which is important 
for light systems, has also been taken into account. There is an excellent 
agreement between the calculations with the modified Glauber model and 
the experimental data. This suggests that the heavy ion reactions in this 
energy range can be explained by the Glauber model in terms of free NN 
cross sections without incorporating any medium modification. 
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I. INTRODUCTION 



The Glauber model (GM) [1] has been employed for describing heavy ion reactions at 
high energies. It is a semiclassical model picturing the nuclei moving in a straight path along 
the collision direction, and gives the nucleus-nucleus interaction [2] in terms of interaction 
between the constituent nucleons (NN cross section) and nuclear density distributions. It 
is a well established model for high energies and has been applied to heavy ion collision 
for describing a number of reaction processes (See e.g. [3-5]). One of the most important 
physical quantities characterizing the nuclear reactions is the total reaction cross section. It 
is very useful for extracting information about the nuclear sizes. The Glauber model has 
been successively used to get the radii of radioactive nuclei from measured total reaction 
cross sections [6]. 

At low energies, the straight line trajectory is modified due to the Coulomb field between 
two nuclei. The Coulomb modified Glauber model (CMGM) [7,8] consists of replacing the 
eikonal trajectory at an impact parameter b, with the eikonal trajectory at the corresponding 
distance Tc of closest approach in the presence of the Coulomb field. Later the non eikonal 
nature of the trajectory around Tc, has also been taken into account [9]. Replacing Tc by 
the distance Ten of closest approach in the presence of both the Coulomb and the nuclear 
field gives the Coulomb plus Nuclear modified Glauber model (CNMGM) [9,10]. This model 
(CMGM/CNMGM) has been widely used in recent literature [11-14]. Warner et. al. [11] 
have shown in their calculations for light projectiles on various targets that trajectory mod- 
ifications are minor. Let us remark that the energies considered by them is much above the 
Coulomb barrier. The Coulomb modified Glauber model has also been used at very high 
energies [12] where trajectory modifications may be ineffective. There have been various 
prescriptions [13,14] to modify the NN cross sections due to nuclear density. Most of the 
work reported earlier [7,8,10,9,13,14] has been done using Gaussian densities. The reaction 
cross section is sensitive to surface density of the colliding nuclei. Depending on how well 
one has produced the surface texture with the Gaussian form, it will make a 5 to 10 % 
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change in the reaction cross section over that done with reahstic densities [15]. 

In the hght of this, we reexamine the various trajectory corrections in the Glauber model 
and calculate the total reaction cross section as a function of energy upto 50 times the 
Coulomb barrier. We quantify the energy range in which a particular correction is effective 
and find that when the center of mass energy of the system becomes 30 times the Coulomb 
barrier, the above modifications to GM are insignificant. Thus, the energy at and above 
which the results of CNMGM and GM coincide depends on the Coulomb barrier between 
the two nuclei and will be different for a light and heavy system. In the present work 
we use exact nuclear densities and free NN cross sections. The center of mass correction 
which is important for light systems has also been taken into account. Comparison of the 
experimental data with our calculations are presented in the plots of aii{E) vs. E. 



Consider the collision of a projectile nucleus A on a target nucleus B. The probabihty 
for occurrence of a nucleon-nucleon collision when the nuclei A and B collide at an impact 
parameter b relative to each other is given by [4,5] 



Here PaO^a) and p^(bB) are the z-integrated densities of projectile and target nuclei 
respectively. t{h)dh is the probability for having a nucleon-nucleon collision within the 
transverse area element dh when one nucleon approaches at an impact parameter b relative 
to another nucleon. All these distribution functions are normalized to one. Here ctjvat is the 
average total nucleon nucleon cross section. 

There can be upto Ax B collisions. The probability of occurrence of n collisions will be 



Here, s = 1 — T{b)aNN- The total probability for the occurrence of at least one NN collision 
in the collision of A and B at an impact parameter b is 



II. THE GLAUBER MODEL 




(1) 




(2) 
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Y:P{n,b) = l-s^^. (3) 



The total reaction cross section an can be written as 

aR^ 277 J bdb (l - s^^) . (4) 

Prom this one can read the scattering matrix as 

\S{b)\'^s^^^{l-T{b)aj,r,)^''. (5) 

In the optical limit, where a nucleon of the projectile undergoes only one collision in the 
target nucleus can be written as 

\S{b)\^c^eM-TmNNAB). (6) 

The scattering matrix can be defined in terms of eikonal phase shift x(6) as 

S(b) = exp (-ix(b)) . (7) 

If we assume ckjvat to be the ratio of real to imaginary part of NN scattering amplitude, the 
eikonal phase shift can be obtained as 

X{b) = ^^NNiaNN + i) AB T{b). (8) 

Here ckjvat will not be directly used for calculating the reaction cross section but will come 
through the correction in the trajectory due to the nuclear field. Once we know the phase 
shift and thus the scattering matrix, we can calculate the reaction cross section and also 
the elastic scattering angular distribution. In the present work we restrict ourselves to the 
calculations of reaction cross section only. 
In momentum space T{b) is derived as 

T{b) = Mqb)SA{cOSBi-ci)fNNiq)qdq. (9) 

Here SA{q) and Ssi—q) are the Fourier transforms of the nuclear densities and Jo{qb) = 
l/27r / exp{—qb cos (f))d(l) is the cylindrical Bessel function of zeroth order. The function 



fNNiq) is the Fourier transform of the profile function t{h) and gives the q dependence of 
NN scattering amphtude. The profile function t{h) for the NN scattering can be taken 
as delta function if the nucleons are point particles. In general it is taken as a Gaussian 
function of width ro as 

exp(-6V2rg) 
'^^^ ^ 2nrl ■ 

Thus, 

/iViv(g)=exp(-roV/2). (11) 

Here, ro is the range parameter and has a weak dependence on energy (see for discussions 
[12]). We use Tq = 0.6 fm [8] in the energy range 2 to 200 MeV. 

The average NN scattering parameter o^jvat is obtained in terms of pp cross section 
(Tpp and nn cross section (Tnm the parameterized forms for which are available in Ref. [8]. 
For projectile energies below 10 MeV/nucleon, we use the prescription given in Ref. [8] for 
computing olnn- For energies above 10 MeV/nucleon, we take the parameterized forms from 
the work of Ref. [20]. These forms produce the elastic scattering data for various systems 
satisfactorily. 

A. Optical potential in the Glauber model 

The usual Glauber optical limit phase shift function is given by 

X{h) j VM{v)dz. (12) 

— oo 

Comparing Eq. (12) with Eq. (8), the optical potential VAf(r) in the momentum representa- 
tion can be identified as, 

AB f 

VN{r) = -fiv ONN{a.NN + i) ^^^^^ j e~''^-''SAiq)SBi-q)fNN{ci)d^q, 

AB r 

= -hv aNNioiNN + i) J jo(qr)SA(q}SB(-q}fNN(ci)q^dq. (13) 

Here jo{qr) is the spherical Bessel function of zeroth order. 



B. Trajectory modifications in the Glauber model 



The basic assumption in the Glauber model is the description of the relative motion of 
the two nuclei in terms of straight line trajectory. For low energy heavy ion reactions, the 
straight line trajectory is assumed at the distance of closest approach Tc calculated under 
the influence of the Coulomb potentials for each impact parameter b as given by, 

r, = (77 + ^ri^ + b^k'')/k, (14) 

which is a solution of the following equation without the nuclear potential VAr(r) 



ZpZre^ b^ 
2n 



E - —- - Rel4r(r) = 0. (15) 



Here 77 = ZpZtg^ /hv is the dimensionless Sommerfeld parameter. The CMGM [7,8] consists 
of replacing the eikonal trajectory at an impact parameter b with the eikonal trajectory at 
the corresponding distance Tc of closest approach. The distance of closest approach Ten in 
the presence of both the Coulomb and nuclear potential [9] is obtained by solving Eq. (15) 
numerically, where Vjv(r) has been derived from Eq. (13). By replacing h by r^n, CNMGM 
is obtained. The non eikonal trajectory [9] around Tc in the presence of the Coulomb field 
is represented by = + (C + l)z^, where the quantity C is given by 

C=^.„. (16) 

Thus in the CMGM, T{b) wiU be simply replaced by T{rc{b))/ ^ {C + I). 

C. The nuclear densities 

The two parameter fermi (2pf) density is given by 

where po — 3/ {At:c^{1 + ^^^)). Thus, the momentum density can be derived [16]as 
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-^U) = ^3 i_e-2. [smx^—-^-xcosxj. (18) 

Where 2; = ndq and x = cq. Here d is the diffuseness and c, the half value radius in terms 
of rms radius R for the 2pf distribution is calculated hy c — (5/3i?^ — T/Sir'^cf — S/Sr^)^/^. 
Here is the proton radius. Equation (9) can be solved numerically for this density and 
the overlap integral can be extracted. In the present work, the density parameters have 
been taken from the compilation of measured charge density distributions [17,18] and are 
tabulated in Table I. 

For lighter nuclei such as ^^C and ^^O, the densities are given in the from of harmonic 
oscillator (HO) densities as 

P{r) = Po (^1 + exp(-^), (19) 

where po = 1/ [1 + l-5a\. The momentum density is given by 

S{q) = po (1 + l-5a - 0.25aq^R'^) exp i -^-^ J . (20) 
Center of mass correction: 

For hghter nuclei such as ^^C and ^^O we also take into account the corrections due to 
Center of motion. Such a correction for harmonic oscillator wave functions is given in Ref. 
[19]. With this the corrected density will become 

S{q) = po (1 + 1.5a - 0.25ag'i?2) exp ( ) exp ( ^ J • (21) 
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TABLES 



TABLE I. Density parameters for various nuclei used in the present work 



A 


Z 


Density 


d/a 


R 






form 


fm 


fm 


12 


6 


HO 


1.082 


1.692 




o 
o 


n\j 




l.Ooo 


28 


14 


2pf 


0.537 


3.150 


40 


20 


2pf 


0.563 


3.510 


58 


28 


2pf 


0.560 


3.823 


90 


40 


2pf 


0.550 


4.274 


208 


82 


2pf 


0.546 


5.513 
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III. RESULTS AND DISCUSSIONS 



The reaction cross sections for light and heavy systems right from near coulomb barrier 
to intermediate energies have been calculated. We calculate the total reaction cross section 
as a function of energy upto 50 times the Coulomb barrier. The data for various systems 
and their references are given in Table II, III and IV. The Coulomb barrier is calculated 
from expression Vc = ZpZtIAA/ {l.b{A^I^ + B^/^)). The reduced mass M is defined as 
M ^ AB/{A + B). 



TABLE II. Reaction cross section data for C on varioTis targets along witli their references 



System 


^labM 




Ref. 




MeV/nucleon 


mb 






9.33 


1444.00 


[21] 


Vc/M = 1.258 


15.00 


1331.00 


[22] 




25.00 


1296.00 






30.00 


1315.00 


[23] 




35.00 


1259.00 


[22] 




15.0 


2165.0 


[22] 


Vc/M = 2.186 


25.0 


2030.0 


!? 




30.0 


2014.0 




12c + 90Zr 


10.0 


2219.0 


[24] 


Vc/M = 3.214 


15.0 


2297.0 


[22] 




25.0 


2415.0 






35.0 


2528.0 




12c + 208pb 


8.0 


1754.0 


[22] 


Vc/M = 5.068 


15.0 


2873.0 






25.0 


3236.0 


55 




35.0 


3561.0 


5? 



9 



TABLE III. Reaction cross section data for °0 on various targets along with their references 



System 


^lab/A 




Ref. 




MeV/nuclcon 


mb 




160 + ^^0 


2.0 


73.86 


[25] 


Vc/M = 1.524 


3.0 


1136.00 






4.0 


1300.00 


5? 




5.0 


1395.00 


11 




9.0625 


1650.00 


[26] 




15.625 


1664.00 


11 




21.875 


1639.00 


11 




30.000 


1655.00 


11 


160 + 28Si 


2.0625 


509.0 


[27] 


Vc/M = 1.90 


2.3750 


765.6 


11 




3.1250 


1341.0 


11 




3.4375 


1451.0 


11 




4.125 


1626.0 


11 




5.0625 


1777.0 


11 




8.9000 


2013.0 


11 




13.45 


2067.0 


11 




94.00 


1757.0 


[10] 
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TABLE IV. Reaction cross section data for °0 on various targets along with their references 



System 






Ref. 




McV/nuclcon 


mb 




160 + ^^Ni 


2.50 


48.2 


[28] 


Vc/M = 2.683 


2.75 


288.7 


5? 




3.00 


539.4 


5? 




3.75 


994.5 


55 




4.375 


1282.0 


55 




5.00 


1485.0 


55 




6.25 


1765.0 


55 




7.50 


1953.0 


55 


16o+208pb 


5.00 


124.1 


[29] 


Vc/M = 5.019 


5.50 


566.6 


55 




6.00 


939.6 


55 




6.50 


1171.0 


[30] 




8.093 


2023.0 


[24] 




12.00 


2847.0 


55 




19.54 


3452.0 


[31] 




94.00 


3600.0 


[32] 
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Figure (1) shows the reaction cross section for ^^O + ^^O system calculated with CMGM 
with and without center of mass (cm.) correction. This correction reduces the reaction 
cross section for lighter systems. Figures (2) to (5) show the reaction cross section for 
^^O on various targets as a function of center of mass energy divided by the Coulomb 
barrier, calculated with the Glauber model (GM), the Coulomb modified Glauber model 
(CMGM) and the Coulomb plus nuclear modified Glauber model (CMNGM) along with 
the data. Figures (5) to (9) show the reaction cross section for ^^C on various targets. 
The center of mass correction has been taken into account in all the calculations. It can 
been seen from all the figures that the most significant effect in this energy range, is by the 
trajectory modification due to the Coulomb field. It is very significant upto energies 6Vc- 
The modification in the trajectory due to the nuclear field tries to bring the results closer 
to GM. We universally find that when the center of mass energy of the system becomes 
30 times the Coulomb barrier, the calculations with the CNMGM match with GM within 
2 to 3%. Thus, the energy at and above which the results of CNMGM and GM coincide 
depends on the Coulomb barrier between the two nuclei and not on the energy per nucleon 
of the projectile. Further, it is observed that there is an excellent agreement between the 
calculations from the modified Glauber model and the experimental data. In contrast to 
[13,14] the present study suggests that the heavy ion reactions in this energy range can be 
explained by the Glauber model in terms of free NN cross sections without incorporating 
medium modification. 

There are higher order corrections [33,34] to the optical limit phase shift function x{b)- 
These corrections are important at higher b but tend to become smaller at larger b [34] 
and are less significant for total reaction cross section which depend mostly upon peripheral 
collisions. They may become important for differential cross sections away from forward 
direction (which probe collisions at smaller b), increasingly at higher energies. The Gaussian 
densities [34] which produce only surface textures of the nuclear density are not adequate 
and realistic densities are to be used (as done in [12]) to calculate these corrections. 
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FIG. 1. Reaction cross section for ^^O + ^^O system calculated with the Coulomb modified 
Glauber model (CMGM) with and without center of mass (cm.) correction. 
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FIG. 2. Reaction cross section for ^^O + ^^O system calculated with the Glauber model (GM), 
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the Coulomb modified Glauber model (CMGM) and the Coulomb plus nuclear modified Glauber 
model (CNMGM) along with the data. 



13 



4000 



Si 




FIG. 4. Same as Fig. 2 but for ^^O + ^^Ni system. 
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FIG. 5. Same as Fig. 2 but for ^^O + ^o^Pb system. 
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FIG. 6. Same as Fig. 2 but for ^^C + ^^C system. 
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FIG. 7. Same as Fig. 2 but for + ^Oca system. 
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FIG. 8. Same as Fig. 2 but for + ^^Zr system. 
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FIG. 9. Same as Fig. 2 but for + ^ospi^ system. 

IV. CONCLUSIONS 

To summarize, we have reexamined the various trajectory corrections in the Glauber 
model and calculated the total reaction cross section as a function of energy upto 50 times 
the Coulomb barrier. The most significant effect in this energy range is by the trajectory 
modification due to the Coulomb field. The modification in the trajectory due to nuclear 
field is also taken into account. We derive V^ir) consistently from Eq. (13). Wc quantify 
the energy range in which a particular correction is effective and find that when the center of 
mass energy of the system is about 30 times the Coulomb barrier, no trajectory modification 
to CM is really required. Exact nuclear densities and free NN cross sections have been used 
in the calculation. The center of mass correction which is important for light systems has also 
been taken into account. There is an excellent agreement between our calculations including 
all the modifications discussed in the manuscript, and the experimental data. In contrast 
to [13,14] the present study suggests that the heavy ion reactions in this energy range can 
be explained by the Glauber model in terms of free NN cross sections without invoking any 
medium modification. It implies that the density effects on the NN cross sections are either 
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absent or very minor. One possible explaination for this is the fact that for heavy ions, the 
contribution to the reaction cross section comes from the surface region where the density 
is very small. 
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